Exact Energy-Time Uncertainty Relation for Arrival Time by Absorption 
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We prove an uncertainty relation for energy and arrival time, where the arrival of a particle at 
a detector is modeled by an absorbing term added to the Hamiltonian. In this well-known scheme 
the probability for the particle's arrival at the counter is identified with the loss of normalization for 
an initial wave packet. Under the sole assumption that the absorbing term vanishes on the initial 
wave function, we show that ATAE > y/ph/2 and {T)AE > 1.37^/ph, where (T) denotes the mean 
arrival time, and p is the probability for the particle to be eventually absorbed. Nearly minimal 
uncertainty can be achieved in a two-level system, and we propose a trapped ion experiment to 
realize this situation. 
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INTRODUCTION 

From Heisenberg's seminal 1927 paper uncertainty 
relations have been recognized as a fundamental feature 
of quantum mechanics. Hcisenberg gives a semi-classical 
heuristic discussion and his "uncertainties" are conceptu- 
ally different in different parts of his paper. Accordingly, 
they are not given a precise quantitative meaning in the 
mathematical language of quantum theory. On the other 
hand, modern textbooks all agree on "the" uncertainty 
relation, namely a version stated and proved by Ken- 
nard 0] in the same year Heisenberg's paper appeared. 
Kennard achieved an important clarification both con- 
ceptually, by defining the uncertainties as the standard 
deviation of operationally well-defined probability distri- 
butions, and also quantitatively, so it becomes possible to 
say that some experiment realizes an uncertainty prod- 
uct within 3% of the absolute minimum. This clarifica- 
tion was so successful that other aspects of Heisenberg's 
paper, like his discussion of the precision of a position 
measurement by microscope versus the momentum dis- 
turbance by the measurement, fell into disrepute. Never- 
theless, these ideas are not only hcuristically meaningful, 
but can be made operationally precise and proved as the- 
orems in the quantum formalism 0, Q • Another group of 
sharp results looks at the position and momentum distri- 
butions like Kennard, but defines the "spread" in a dif- 
ferent way. Particularly interesting is a non-parametric 
version in terms of entropies [5|— |TJ , which has applica- 
tions in cryptography Q. Since many experiments are 
nowadays approaching quantum limits, it is perhaps be- 
coming more important to think of uncertainty relations 
not as a single fact, but as a circle of ideas in which tradi- 
tional heuristic interpretations coexist with an increasing 
family of exact results. 

The ambiguities in the literature multiply when we go 
from position-momentum uncertainty to energy-time un- 
certainty. In most cases, energy-time uncertainty rela- 
tions are invoked in a handwaving fashion only. There 



are very few conceptually clear and quantitatively mean- 
ingful formulations (see [H for a review). For one kind 
of arrival time observables, the so-called covariant ob- 
servables (lp| - [l2l |. a Kennard-like uncertainty relation is 
well known: Then AE is the standard deviation of the 
energy observable and AT the standard deviation of ar- 
rival time. However, this notion of arrival times is rather 
inflexible and makes no sense for finite dimensional sys- 
tems. These drawbacks can be cured in a more realistic 
theory of arrival times [l3|, 14 1, in which the counter is 



modeled by an absorbing term in the Hamiltonian. 

The purpose of this paper is to provide a sharp quanti- 
tative formulation of uncertainty for energy and absorp- 
tive arrival time. The possibility of such an uncertainty 
relation is somewhat surp rising (and totally missed in the 



literature, including [14|), since the dynamics including 



the counters, which is used to define AT, is not generated 
by the Hamiltonian used to define AE, so there would 
appear to be no universal trade-off inequality. Neverthe- 
less, we will prove such a relation under only a mild and 
natural assumption on the initial state. An immediate 
advantage of the absorptive arrival time approach is that 
it also applies to certain finite dimensional models used 
in quantum optics. Thus, we propose a concrete minimal 
uncertainty experiment with trapped ions. 

Our paper is organized as follows. After a brief in- 
troduction to absorptive arrival times, we will state the 
relation and the conditions for near-equality. The main 
ideas of the proof are sketched in the following section, 
with a full proof given in the Appendix. We then show 
that (very nearly) minimum uncertainty can be realized 
in a standard quantum optical setting. 



ARRIVAL TIMES 

Consider a quantum system with Hilbert space Ji and 
Hamiltonian H . Starting from some initial state ip £ "H, 
we would like to determine the probability distribution 
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of arrival times at some counter. There are different ap- 
proaches to this problem, varying in the degree of detail 
with which the counter is described. 

The coarsest, and simplest, description focuses just on 
transformation behavior: starting from the time-evolved 
state ipt = exp(—iHt/h)ip we should get the same ar- 
rival time distribution, but shifted by t. Covariant ar- 
rival observables with this property have been studied 
extensively 0,1 121 llEfl. and satisfy a general energy-time 
uncertainty relation However, transformation be- 

havior alone is not sufficient to single out a convincing 
model for a given experimental counter array. Moreover, 
this approach requires the Hamiltonian to have a purely 
continuous spectrum and is hence limited to infinite di- 
mensional Hilbcrt spaces. 

At the other extreme we can make a detailed detec- 
tor model, e.g., by interaction with an ionizable atom 
The drawback of this scheme is that, although un- 
doubtedly physically correct, the interacting system is 
very hard to treat, and a disproportionate amount of the 
analysis of a given experiment would go just into the de- 
tection process. 

The absorptive arrival times approach adopted in this 
paper is a good compromise between these extremes: 
It was first worked out in detail in [l3| and describes 
the detector by a non-hermitian term — iD added to 
the Hamiltonian H, which is thus replaced by K = 
H — iD. Thereby the unitary time evolution operator 
Ut = exp(—iHt/H) is modified to a semigroup of contrac- 
tions, i.e., operators B t = cxp(—iKt/h) (t > 0) such that 
1 1 St 1 1 < 1. We now interpret the loss of normalization, 
i.e., 1 — H-BtV'l! 2 , °f a quantum state ip as the probability 
that the particle did arrive before time t. More generally, 
for any time interval [t, s] with < t < s the probability 
of arrival in that interval is given by the expectation of 
the operator 



F([t,. 



BTB t - B*B S 



(1) 



Thus, given a quantum state ip, we get a probability den- 
sity P(f) on the positive time axis by 



P(t) = -- ^(BrtW), 
p at 

p = 1 — lim ||-BtV>|| 

f— >oo 



(2) 
(3) 



Here the limit R = linit-^oo B^B t exists because it is over 
a decreasing family of positive operators, and p is the to- 
tal absorbtion probability. For real detectors not every 
absorption actually leads to a detected click. Suppose the 
probability for this process is q. Then we would observe 
a click in [t, t + dt] with probability qpP(t)dt. Hence F(t) 
can be determined from the experimental data by nor- 
malizing the observed click distribution, independently 
of q or p. The important distinction between these two 
probabilities is that the process leading from absorption 
to detection is independent of the particle dynamics and 



introduces no change in the further quantum evolution of 
the particle. Therefore, q cannot enter in the uncertainty 
relation, but as we will see, p does. 

We denote by (T) and (T 2 ) the first and second mo- 
ment of the probability distribution P(£) after and 
set (AT) 2 = (T 2 ) - (T) 2 . Defining (AE) 2 = (ip\H 2 iP) - 
(iplHip) 2 , we can hence look for a universal lower bound 
on the product AT • AE. 

Without further conditions such a lower bound cannot 
hold. Indeed, AT can be computed knowing K and ip, 
whereas AE depends on H and ip. For example, if we 
now set K = H — ictl, we get p = 1 and P(i) = 2ae~ 2at 
independently of ip. Clearly, this cannot imply any con- 
straint on the energy distribution. Similarly, if the initial 
state is a joint eigenvector of H and K, then AE = 
and AT has a finite value belonging to an exponential 
distribution. 



UNCERTAINTY RELATION 

From the last paragraph it seems that it makes no sense 
to look for a general time-energy uncertainty relation in 
this setting. However, as we will now show, a simple 
and physically natural condition suffices to derive one. 
Loosely speaking, the condition is that the initial wave 
function has no overlap with the detector. More formally, 
if D describes the detector as explained above, we want 
Dip = 0, or Hip = Kip. Since these operators are usually 
unbounded, we also have to specify the domains. Writing 
dom X for the domain of the operator X , we require that 



tp £ dom K n dom H and Hip = Kip. 



(4) 



The main result of our paper is that under this condition 



AT • AE > - y/p 



(5) 



The dependence on p expresses the fact that for a small 
detection operator D only a few particles are ever de- 
tected (p»0), so observing arrival times cannot imply a 
strong constraint on AE. The power v /p is explained in 
the sketch of proof below. 

The arrival time distribution P(i) is always supported 
by the positive time axis R+. Therefore, the mean arrival 
time (T) is always positive, and can take the place of 
AT in the uncertainty relation. Often (T) is of more im- 
mediate relevance than AT. Consider, for example, the 
decay of a metastable state. The initial wave function 
is trapped inside a potential barrier, through which it 
will eventually tunnel. A detector is placed at a distance 
from the potential (so our initial condition Kip = Hip is 
satisfied). Then (T) is directly an expression of the life- 
time of the metastable state. Under the same conditions 
as for ([5]) we prove that 



(T) • AE > C h y/p 



(6) 
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where C = 2(— Zi/3)^ 3 / 2 ' « 1.376 is a numerical constant 
involving the first negative zero of the Airy function. 

Cases of (nearly) minimal uncertainty are in both cases 
connected to specific probability distributions P m in(i)- 
Let (1 + e) denote the ratio of the left hand side to the 
right hand side in ([5} or ([6|) . Then, for a suitable scaling 
factor A and shift r. 



dt P(i) - AP mi „(A(i - r)) < 



(7) 



For ([5]), Pmin(^) is Gaussian, and 7 = -\/2, as for the stan- 
dard position- momentum uncertainty relation. For ([5]), 
Pmin(i) is the square of the Airy function, and 7 = 1.888. 
We remark that no such conclusion can be drawn for the 
energy distribution. In fact, although ([5]) is nearly sharp 
for the two-level system discussed below, the energy dis- 
tribution is concentrated on two points, and is hence very 
different from a Gaussian. Another consequence is that 
equality (i.e., e = 0) is impossible in ([5]), because every 
Gaussian has a tail extending to the negative half axis. 
Equality is possible for however. 



SKETCH OF PROOF 

The basic idea of the proof is a so-called dilation con- 
struction, by which the system is transformed in such a 
way that energy and time become conjugate self-adjoint 
operators, and theorems about position-momentum pairs 
can be used. The full argument allowing also unbounded 
H and K is given in the Appendix. Here we give a rough 
version in theoretical physics (rather than mathematical) 
style. For simplicity we focus on the AE AT relation (|5|), 
and assume that K = H — iD for a bounded positive op- 
erator D. For the dependence on the total absorption 
probability p we only give an argument showing why the 
square root is the expected power. For the main part of 
the proof we then just assume that p = 1 for all initial 
states. This is true, e.g., in typical finite dimensional 
systems. 

To understand the ^Jp factor in ([5]) and (|6]), consider 
some T-L, H, K,tp as above, with full absorption p = 1. 
Compare it with the following modification, which just 
adds a part of the Hilbert space that is "not seen" by 
the detector. Formally, the Hilbert space T~L is enlarged 
by an additional orthogonal summand, on which K van- 
ishes, and which contains an eigenvector <pQ of H with 
(4>o\H<fio) — (ip\Hip) = E. In this system we con- 
sider the initial vector ip' = y/p ip + s/1 — p' </>o , where 
< p' < 1. This p' is precisely the absorption parame- 
ter for the extended system. The time distribution does 
not change at all since we defined it as normalized by 
p' . The energy expectation also does not change, but 
A^/E 2 = p'A^E 2 + (1 - p'XA^E) 2 = p'A^E 2 . Hence 
([S]) and (O, which hold with p = 1 for ip, are equivalent 



to their versions for ip' with the factor y/p 1 . More com- 
plex ways in which the particle might not be seen by the 
detector are covered by the full proof. From now on we 
just assume p = 1. 

We will associate with any wave function ip € H an- 
other wave function ip, which is a function of time, so that 
|V>(£)| 2 is the arrival probability density. In other words, 
ip is a wave function in a time representation. For each 
t, ip(t) lies in the original Hilbert space H. We define 
ip = Jip with the linear operator 



(JiP){t) 



/ 2/HD 1 / 2 B t iP if t > 
otherwise 



(8) 



where the square root of the positive operator D is taken 
in the operator sense. Then, for t > 0, 

\\{Jm)f = \{B t ip\DB t iP) 

= -j t {B t iP\B t iP)=¥{t) . (9) 

Here we used the assumption that p = 1 for all states, so 
the normalizing factor p in ([2]) can be omitted. Equiva- 
lently, J* J = I. We denote by T the self-adjoint oper- 
ator of multiplication by t. Then the required moments 
ofP(i) are (T) = $\fip) and (T 2 ) = (^|f 2 ^). 

The translation in the time representation is generated 
by the operator H = ihdj (dt). By differentiating ((5J and 
using the assumption Hip = Kip, we find, for t > 0, 



(HJtp)(t) = ihy/2/HD 1/2 B t {-iK/h)ip 
= (JA»(<) = (Jfty)(t). 



(10) 



This obviously also holds for t < 0. At t = the def- 
inition ([5]) would appear to allow a jump discontinuity, 
leading to a 8(t) contribution in (fTOj) . But since Dip = 
we also have (Jip)(0) = 0, and hence there is no jump. 
We can thus take (fT0|) as an equation for functions, i.e., 
Hip = HJtp = J Hip. Therefore, using J* J = 1, the vari- 
ances AE 2 = \\Hif>\\ 2 — (ip\Hip) 2 and the corresponding 
AE 2 = \\Hip\\ 2 — (ifj\Hip) 2 are the same. The time-energy 
uncertainty relation ([5]) thus follows from the standard 
position-momentum one, applied to "position" T and 
"momentum" H. 

It is well-known that the uncertainty relation can be 
derived from the operator inequality H 2 + T 2 > 1 and a 
dimensional scaling relation. A similar argument applied 
to wave functions only on the positive time axis, using 
the inequality H 2 + T > XI gives ([6]). A in this equa- 
tion is determined from the ground state problem of a 
particle in a linear potential in front of a wall, which is 
solved in terms of the Airy function. Using the gap to 
the first excited state in this problem (and the oscillator 
Hamiltonian) gives the estimates for near-equality ([7]). 
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P(t), Pmin(t) 




FIG. 1. Thick curve: Arrival time probability density P(i) 
after for fl/2 = h = 1 and 7 = 2v / 2- Dashed: probability 
distribution P m i n (i). 



EMISSION FROM A TWO-LEVEL SYSTEM 

Let us consider the simplest possible system to which 
the relations (|5|) and ([6]) apply. Unlike the covariant ob- 
servable approach, in which the spectrum of H has to be 
continuous, finite dimensional systems are included. We 
consider a two-level system with 



H = 



2 \n 



D = 



fi/o 



2 10 7 



(11) 



where fl, 7 > are parameters. The relevant quantity is 
7/0, since we could make 0/2 = H = 1 by a choice of 
units. As the initial state we take tp = (1,0) T , so our 



iD)ijj = Kijj 



so 



basic assumption Hip = (0, M7/2) T = (H — 
is satisfied. Obviously, (H) = and (H 2 ) 

ae = fa 

One can explicitly exponentiate B t = exp(—iKt/h), 
and hence compute the probability density P(i). The 
result is shown in Fig. [TJ The moments are also readily 
calculated. Both (T) and AT attain their minimum when 
7 = \/2Q. Then the uncertainty inequalities are satisfied 
as 1/V2 w 0.707 > 0.500 for ©, and a/2 « 1.414 > 1.376 
for (O, so (T)AS reaches the minimum to within 3%! 
This implies by ([7]) that the arrival time distribution P(t) 
must be close to P m i n (i)- The comparison is shown in 
Fig.ffl 



EXPERIMENTAL IMPLEMENTATION 

Instances of the uncertainty inequality as such are not 
an interesting experimental target, since it is impossible 
not to implement it. On the other hand, instances with 
near minimal uncertainty can be a challenge. The two- 
level scheme approximations described above naturally 
arises in quantum optics when looking at an atom in- 
teracting with lasers. A feasible level scheme is the one 
of a single 40 Ca + atom (similarly level schemes exist in 
other atoms), trapped in a Paul trap, which is shown in 



Fig. [2j The transitions 1—^2 and 2 — s- 3 are driven 
by on-resonance lasers resulting in Rabi-frequencies O12 
resp. O23. The narrow quadrupole transition 1 — > 2 
can be frequency resolved by a narrow-linewidth laser, 
whereas the 2 — > 3 transition is selected by choosing a- 
polarized light. If | O23 1 -C r3_>4, we will get approxi- 
mately an efficient decay rate 7 = |023| 2 /r3->4 for level 
2. If we consider only levels 1 and 2, the corresponding 
effective Hamiltonian, within the rotating wave approx- 
imation and in the interaction picture to get rid of any 
time dependence, will be K = H — iD as in (fTTI) . with 
this 7 and ft = ili2- The initial level 1 can be efficiently 
prepared via optical pumping. The arrival time in state 
2 is then measured by measuring the time of the first 
spontaneously emitted photon on the 3 — > 4 transition. 
In the proposed scheme, the total absorbtion probability 
p is smaller than 1 due to spontaneous scattering events 
from state 3 back to state 2. However, these events are 
suppressed by more than a factor of 200 due to favorable 
branching ratios and small Clcbsch-Gordan coefficients 
on the 2—^3 transition. All other spontaneous decay 
channels result merely to a reduced detection efficiency 
and do not affect the uncertainty relation. A typical Rabi 
frequency of O12 ~ 27r x 100 kHz would require an easily 
achievable ^23 ~ 27r x 1.73 MHz to reach approximately 
the minimum of the product (T)AE, see Fig. [TJ The 
required time resolution for the photon detector in this 
scenario is on the order of a few ns, well within reach 
of current technology. Alternative implementations in- 
clude vacuum-stimulated Raman transitions in an atom 
strongly coupled to a leaky cavity 



17] 




FIG. 2. Level scheme of Ca , the relevant levels 1,2 for 
a near-minimum uncertainty experiment are emphasized, the 
transitions with their branching ratios j3 are: Ai^2 = 729.147 
nm, F 2 ^i = 2tt x 0.0544 Hz, /3 2 ->i = 1; A 2 ^s = 854.209 nm, 



r 3 ^ 2 = 2tt x 105 kHz, 
T 3 ^4 = 2vr x 21.2 MHz, 



#3- 
03- 



■ 2 = 0.87o: 
> 4 = 92%. 



A4- 



393.366 nm, 
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OUTLOOK 

The approach given here can be applied to other situa- 
tions involving time in quantum mechanics, such as tun- 
neling and decay scenarios. A standard detector model 
for such applications is the pointlike counter, formally 
given by D = XS(x). In the finite dimensional (n > 2) 
setting it remains an interesting problem to characterize 
the arrival time densities P(f), which can be engineered 
by a suitable choice of laser couplings and initial states. 
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APPENDIX: PROOF 

The formal proof of the statements in the paper fol- 
lows the sketch given there, with the following main dif- 
ferences: 

• The operators A and H are allowed to be un- 
bounded, and domain questions are treated explic- 
itly. This covers cases, where A — H is not mean- 
ingful as an operator (e.g., a 5-function potential). 
Therefore, the difference D = i(K — H) appears 
nowhere, and the dilation operator J has to be de- 
fined differently from ((SJ. 

• There is no assumption of total absorption, i.e., 
p < 1 is allowed. 

• The differentiation (|I0[) is only carried out under 
suitable scalar products, avoiding the discussion 
distributional derivatives. 

• We set h = 1 for notational convenience. 



The Dilation 

We begin by recapitulating the dilation construc- 
tion from Since B t is a contraction semigroup, 
—d/dt(B t 4 l \Bt4)) is a positive semidcfinitc scsquilincar 
form on dom A. We denote its completion by £, which 
we call the exit space, and by j : dom/v — > £ the embed- 
ding map, so that for 0, tp S dom A: 



(m4>)e = -j t (B t 4>\Bt 
= i{{i>\K<t>) 



The basic dilation operator J maps into the space 

H = C 2 (R,dt;£), (13) 

which is the space of "^-valued wave functions" t H> 
*(t) e £ such that ||4<|| 2 = /~<ft||*(t)||| < oo. For 
tp G dom.fr, we define 



(JV#) 



Then 



\W\\ 



j{B t il>) if £ > 
if t < " 



dt {jB t tj)\jBti)) 



f°° ( d 



t=0 

{i/>\(l-R)i/>), 



(14) 



where R — lim B^B t . 

t— >oo 



(15) 
(16) 



Hence \\Jip\\ < || "011 1 an( i J extends by continuity to a 
unique operator J : % —> C 2 (R, dt; £) , which we will 
denote by the same letter. 



Time Distributions 

Now in the dilation space H we have a time operator T, 
which acts by multiplication with t. Its spectral projec- 
tions F([t, s]) are a dilation of the arrival time observable 
F. That is, in analogy with (| 15[) we get, for < t < s 
and cj), if) € dom A': 

(Jtf>\F([t,s])Jip) = [ dr (j B T <f>\j B T i/)) 



= -(B r ct>\B r il>) 



\F([t,sM . 



(17) 



Since the operators involved are bounded and dom A is 
dense, this is summarized in the operator equation 



F([t,. 



J*F([t,s])J. 



(18) 



=0+ 



(12) 



In particular, the probability density P after ([2|) is equal 
to P(i) = \\(Jip)(t)\\ 2 /p, i.e., the probability density as- 
sociated to the time operator T by the vector 



V> = — Ji>. (19) 
This vector is normalized, because J* J = 1 — R, so p = 

WW). 
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Energy Distributions 

The dilated energy operator H is defined as the gener- 
ator of the translations Ut = exp(-iHt), where 



(20) 



for t,r £ R. Obviously, for t,T > and tp £ domK, 
we have (U T Jip)(t) = j(B t + T ip) = (JB T ip)(t), and hence 
J*UtJ = J*JB t . By taking adjoints we also get an ex- 
pression for t < 0. To summarize, 



J*U t J 



(l-R)B t if t > 
B* t (l-i?) if i < 



(21) 



Again, these expressions are directly verified on domK 
and extended by continuity to all of H . This shows that 
in the special case R = 0, U is a dilation of the semigroup 
Bt in the sense of Sz.-Nagy [2~~ 



21] 



Since B* t RB t 



R, 



we can commute generators and R in the sense that 

R dom K C dom K* , 
RK = K*R . 



(22) 



The ^-distribution with respect to ip from (|19p is the 
measure (ip\E(du))ip), where E is the spectral measure 
of H . The characteristic function (Fourier transform) of 
this distribution is 



C(t) 



P 



e- iut {^j\E{dw)^) = (iP\U t iP) 

O 

(1>\J*U t Jil>), 



(23) 



which can be evaluated further using (|21[) . The moments 
of the iJ-distribution are obtained by differentiating C 
at zero. In general, of course, the piecewise defined func- 
tion (f2~Tj) is not differentiable at zero. However, with our 
assumption on ip, we can establish that the right and left 
derivatives coincide. 

The conditions ip £ dom K n dom H and Kip = Hip 
imply that \\j(i>)\\ 2 = i{{i'\Hi')-{Hip\ip)) = 0, so jip = 
for such i/j. Moreover, for all <fi £ domK, 

WW) - (m<p) = -iw>\j<j>) = o, 

whence K*tp = Kijj = Hip. Now the first derivative of 
the characteristic function becomes 

tpC ({J > - { (iP\K*(l - R)ij>) for t = 0- {M > 
= (ip\(l - R)Hijj) = (Hijj\(l - R)iP). 

Similarly, for the second order, and with the additional 
property ip £ dom_fC 2 we get 

-«r'Vfrt-J Wf 1 -^) fort = 0+ 
1 ' ~ I (^\(K*) 2 (1- R)ip) fort = 0- 



Hcnce, with A = (ip\Hip) we get the variance 

(AE) 2 = ((H - AIty| (IT - M)ip) 

> ((H - Xl)i'\(l - R)(H - Xl)iP) 
= p(-C"(0) - 2i\C'(0) + X 2 C(0)) 

= p I (w- A) 2 tr (ip\E(duj)ip) 



> pmm J (u - A) tr{ip\E(du)ip) 



To summarize: 



AE > y/p AH. 



(26) 



(HiP\(l- R)HiP) 



(25) 



In particular, AEAT > ^/pAHAT > ^/p/2, which 
proves ([5]). However, since we also want (|6]) and the error 
estimate (|7jl . we need a more detailed analysis. 



Conversion to ground state problems 

The inequalities ([5]) and ^ are for a product of mo- 
ments. It is advantageous to turn this into inequalities 
for the expectations of a sum of operators. For this we 
use a standard trick, which basically amounts to a di- 
mensional analysis, and hence assures that we estimate 
quantities of the physical dimension time x energy. If we 
change the unit of time, we transform T i— > rj~ 1 T, and 
H i— > r/H. This change is also implemented by a unitary 
operator in R. Similarly, we can shift H and T by mul- 
tiples of the identity with unitary operators. This is the 
idea behind choosing the operators 

X = rf{H -e) 2 + 7f 2 (f -t) 2 , (27) 
Y = r, 2 {H - ef + r i - 1 f, (28) 

where Y is considered as a quadratic form on the sub- 
space of functions vanishing for t < 0. These are just 
standard Schrodinger operators with H interpreted as 
the momentum and T as the position. In this guise they 
are well-known: X is the Hamiltonian of the harmonic 
oscillator, and Y describes a particle in a linear poten- 
tial in front of a wall. Here the quadratic form point 
of view is dictated by our aim to bound the expectation 
value (ip\Yip) for functions in the positive-time subspace, 
and corresponds in operator theoretic terms to taking 
the Dirichlet boundary condition at t = 0. Note that 
this also corresponds to our discussion preceding ([24| . 
which says that for the initial vectors ip satisfying our 
assumptions, we have jip = ip(0) = 0. Both operators 
have purely discrete spectrum, because the potentials di- 
verge at infinity. Moreover, the eigenvalues do not de- 
pend on the parameters r, e, 77, because of the unitary 
equivalence. To compute the bottom eigenvalue, we can 
thus set e = 0, r = 0, 77 = 1. 
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For X we have the eigenvalues x n = (2n+l) (note that 
we left out a conventional factor of 1/2 from the definition 
of the oscillator Hamiltonian). The ground state is, of 
course, a Gaussian. 

The eigenvalue equation for Y becomes the differential 
equation 

-r(t)+m) = vm, (29) 

which is solved in terms of the Airy function Ai by 



(30) 



Indeed, Ai is defined to be that solution of — Ai"(t) + 
fAi(f) = 0, which is square integrable at t — > +00. The 
eigenvalue equation comes out of the boundary condition: 
we must have ip(0) = Ai(0 — y) = 0, i.e., (— y) must be 
a zero of Ai on the negative half axis. These zeros are 
well-known and tabulated [22|, Table 9.9.1]. The first two 
are 



Z x = -y = -2.33810 74104, 
Z 2 = -yx = -4.08794 94441. 



(31) 
(32) 



To get inequality ([5]) we choose e = (^}\Hi/j) and r = 
(T). Taking expectations of the inequality X > xol we 
find 



rfAH 2 +ij- 2 AT 2 > x = l. 



(33) 



Minimizing this expression over 77, i.e., setting r/ 2 = 
AT /AH we get 2 AT AH > 1 for the canonical pair 
T, iJ, and (O as stated by using (f2l)]) . We went through 
this well-known argument just to stress the analogy with 
the proof of (|6|), for which the corresponding expression 
using Y reads 



V 2 AH 2 + V - L (T) >y . 



(34) 



This time the minimum is attained at rf = (T) / (2AE 2 ), 
giving 



/ (T)AH 2 \ 3 
3 P-^ ^Vo- 



(35) 



Together with (|2"6"]l this gives ©■ It now remains to 
show ([7]) for the almost minimal case. For this we have 
to first provide an elementary estimate, which quantita- 
tively captures the heuristic idea that in a gapped ground 
state problem only states close to the ground state can 
have expectations close to the bottom eigenvalue. 



Lemma on gapped ground states 

Let A be an operator, which has a lowest, non- 
degenerate eigenvalue ao with normalized cigenfunction 



4>o, so that the rest of the spectrum which lies in the 
half-axis [04,00) with gap ax — ao > 0. That is, we can 
write 



A = a \(j) ) (<f) \ + Ax with 
Ai>ai(I-|0oXfo|). 
Now suppose that for some state p we have 

tr pA < ao + Oi, 

where a > 0. Then we claim that 



|p-|0o)<0o|||i<2 



ax - ao 



(36) 
(37) 



(38) 



(39) 



where ||A||i = tr %/ X*X denotes the trace norm of the 
operator X. The reason to write this in terms of density 
operators and trace norms (although we only need the 
case of a pure state p = later on) is that in this 

form it is clearer that the bound is equivalent to 



\tx{ P S) - (MSfo)\ < 2 



ax - a 



(40) 



for all operators S with 11511 < 1. Note that for any 
two density operators we have \\pi — P2W < 2. Hence 
the bound trivializes for a > ax — ao, but gives some 
information as long as the expectation ao + a lies in the 
gap. 

For the proof we evaluate the inequality 
A > ail - (ax - a Q )\(/)o) ((f>a\ 
with p, getting 

ao + a > txpA > ax - (ax - a o )(^o|p0o), 

which amounts to an estimate for the fidelity (0o|p0o) 5 
namely 



l-(0oW>o) < 



ax - a 



(41) 



The trace norm estimate follows from this by (23l Lemma 
4.1]. 



Almost minimal case 

Let us begin with the variance inequality. Of course, 
it is a standard exercise in almost every quantum physics 
textbook to show that the uncertainty relation holds with 
equality exactly for Gaussians. Much less well-known is 
the stability statement for this minimal case, which we 
prove here. It makes an almost equally simple exercise, 
so one might have expected to see it posed in quantum 
mechanics courses the world over. However, we could not 
find a reference for it. 
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So suppose that ATAE < (1 + e)y/p/2. By J2BJ this 
implies 2 AT AH < (1 + e). The left hand side of this 
inequality is just the expectation of X with the minimiz- 
ing choices of constants. Hence the assumption (|55|) is 
satisfied and the estimate (|4TI|) applies, which gives 



|tr(pS)-<0o|Sto)|<2 



2 5 



(42) 



for any S* with < 1. We choose S as the multiplica- 
tion operator with s(t) = sign(P(i) — P m i n (t)). Then 



-mi„||l = / dt |P(t)-Pmin(t)| 
/0 



dt s(i)(P(*)-PmmW) 



< a/27. 



(43) 



The argument for near equality of ([6]) is completely 
analogous. Suppose that (T}AE < (1 + e)c^/p, so 
(f)AH < (1 + e)c. Then the expectation of F, which 
is the left hand side of ([35]) is by a factor (1 + e)( 2 / 3 ) 
larger than the absolute minimum. We can hence apply 
the estimate with a = ((l+e)( 2/3 )-l)z/ < 2y e/3, which 
gives 

(44) 

3(yi - yo) ^ 
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